Building up a genomic surveillance platform for SARS-CoV-2 in the middle of a pandemic: a true North-South collaboration.

Next-generation sequencing technology has revolutionised pathogen surveillance over the last two decades. However, the benefits are not equitably distributed, with developing countries lagging far behind in acquiring the required technology and analytical capacity. Recent declines in the cost associated with sequencing-equipment and running consumables have created an opportunity for broader adoption. During the COVID-19 pandemic, rapid diagnostics development and DNA sequencing revolutionised the ability to diagnose and sequence SARS-CoV-2 rapidly. Socioeconomic inequalities substantially impact the ability to sequence SARS-CoV-2 strains and undermine a developing country's pandemic preparedness. Low- and middle-income countries face additional challenges in establishing, maintaining and expanding genomic surveillance. We present our experience of establishing a genomic surveillance system at the Aga Khan University, Karachi, Pakistan. Despite being at a leading health sciences research institute in the country, we encountered significant challenges. These were related to collecting standardised contextual data for SARS-CoV-2 samples, procuring sequencing reagents and consumables, and challenges with library preparation, sequencing and submission of high-quality SARS-CoV-2 genomes. Several technical roadblocks ensued during the implementation of the genomic surveillance framework, which were resolved in collaboration with our partners. High-quality genome sequences were then deposited on open-access platforms per the best practices. Subsequently, these efforts culminated in deploying Pakistan's first SARS-CoV-2 phyllo surveillance map as a Nextstrain build. Our experience offers lessons for the successful development of Genomic Surveillance Infrastructure in resource-limited settings struck by a pandemic.

Luna Virus and Helminths in Wild Mastomys natalensis in Two Contrasting Habitats in Zambia: Risk Factors and Evidence of Virus Dissemination in Semen.

Transmission dynamics and the maintenance of mammarenaviruses in nature are poorly understood. Using metagenomic next-generation sequencing (mNGS) and RT-PCR, we investigated the presence of mammarenaviruses and co-infecting helminths in various tissues of 182 Mastomys natalensis rodents and 68 other small mammals in riverine and non-riverine habitats in Zambia. The Luna virus (LUAV) genome was the only mammarenavirus detected (7.7%; 14/182) from M. natalensis. Only one rodent from the non-riverine habitat was positive, while all six foetuses from one pregnant rodent carried LUAV. LUAV-specific mNGS reads were 24-fold higher in semen than in other tissues from males. Phylogenetically, the viruses were closely related to each other within the LUAV clade. Helminth infections were found in 11.5% (21/182) of M. natalensis. LUAV-helminth co-infections were observed in 50% (7/14) of virus-positive rodents. Juvenility (OR = 9.4; p = 0.018; 95% CI: 1.47-59.84), nematodes (OR = 15.5; p = 0.001; 95% CI: 3.11-76.70), cestodes (OR = 10.8; p = 0.025; 95% CI: 1.35-86.77), and being male (OR = 4.6; p = 0.036; 95% CI: 1.10-18.90) were associated with increased odds of LUAV RNA detection. The role of possible sexual and/or congenital transmission in the epidemiology of LUAV infections in rodents requires further study, along with the implications of possible helminth co-infection.

Metagenomic Pathogen Sequencing in Resource-Scarce Settings: Lessons Learned and the Road Ahead.

Metagenomic next-generation sequencing (mNGS) is the process of sequencing all genetic material in a biological sample. The technique is growing in popularity with myriad applications including outbreak investigation, biosurveillance, and pathogen detection in clinical samples. However, mNGS programs are costly to build and maintain, and additional obstacles faced by low- and middle-income countries (LMICs) may further widen global inequities in mNGS capacity. Over the past two decades, several important infectious disease outbreaks have highlighted the importance of establishing widespread sequencing capacity to support rapid disease detection and containment at the source. Using lessons learned from the COVID-19 pandemic, LMICs can leverage current momentum to design and build sustainable mNGS programs, which would form part of a global surveillance network crucial to the elimination of infectious diseases.

Workshop-based learning and networking: a scalable model for research capacity strengthening in low- and middle-income countries.

Science education and research have the potential to drive profound change in low- and middle-income countries (LMICs) through encouraging innovation, attracting industry, and creating job opportunities. However, in LMICs, research capacity is often limited, and acquisition of funding and access to state-of-the-art technologies is challenging. The Alliance for Global Health and Science (the Alliance) was founded as a partnership between the University of California, Berkeley (USA) and Makerere University (Uganda), with the goal of strengthening Makerere University's capacity for bioscience research. The flagship program of the Alliance partnership is the MU/UCB Biosciences Training Program, an in-country, hands-on workshop model that trains a large number of students from Makerere University in infectious disease and molecular biology research. This approach nucleates training of larger and more diverse groups of students, development of mentoring and bi-directional research partnerships, and support of the local economy. Here, we describe the project, its conception, implementation, challenges, and outcomes of bioscience research workshops. We aim to provide a blueprint for workshop implementation, and create a valuable resource for bioscience research capacity strengthening in LMICs.

Chemogenomics identifies acetyl-coenzyme A synthetase as a target for malaria treatment and prevention.

We identify the Plasmodium falciparum acetyl-coenzyme A synthetase (PfAcAS) as a druggable target, using genetic and chemical validation. In vitro evolution of resistance with two antiplasmodial drug-like compounds (MMV019721 and MMV084978) selects for mutations in PfAcAS. Metabolic profiling of compound-treated parasites reveals changes in acetyl-CoA levels for both compounds. Genome editing confirms that mutations in PfAcAS are sufficient to confer resistance. Knockdown studies demonstrate that PfAcAS is essential for asexual growth, and partial knockdown induces hypersensitivity to both compounds. In vitro biochemical assays using recombinantly expressed PfAcAS validates that MMV019721 and MMV084978 directly inhibit the enzyme by preventing CoA and acetate binding, respectively. Immunolocalization studies reveal that PfAcAS is primarily localized to the nucleus. Functional studies demonstrate inhibition of histone acetylation in compound-treated wild-type, but not in resistant parasites. Our findings identify and validate PfAcAS as an essential, druggable target involved in the epigenetic regulation of gene expression.

The Plasmodium falciparum ABC transporter ABCI3 confers parasite strain-dependent pleiotropic antimalarial drug resistance.

Widespread Plasmodium falciparum resistance to first-line antimalarials underscores the vital need to develop compounds with novel modes of action and identify new druggable targets. Here, we profile five compounds that potently inhibit P. falciparum asexual blood stages. Resistance selection studies with three carboxamide-containing compounds, confirmed by gene editing and conditional knockdowns, identify point mutations in the parasite transporter ABCI3 as the primary mediator of resistance. Selection studies with imidazopyridine or quinoline-carboxamide compounds also yield changes in ABCI3, this time through gene amplification. Imidazopyridine mode of action is attributed to inhibition of heme detoxification, as evidenced by cellular accumulation and heme fractionation assays. For the copy-number variation-selecting imidazopyridine and quinoline-carboxamide compounds, we find that resistance, manifesting as a biphasic concentration-response curve, can independently be mediated by mutations in the chloroquine resistance transporter PfCRT. These studies reveal the interconnectedness of P. falciparum transporters in overcoming drug pressure in different parasite strains.

Estimation of Secondary Household Attack Rates for Emergent Spike L452R Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variants Detected by Genomic Surveillance at a Community-Based Testing Site in San Francisco.

BACKGROUND: Sequencing of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral genome from patient samples is an important epidemiological tool for monitoring and responding to the pandemic, including the emergence of new mutations in specific communities. METHODS: SARS-CoV-2 genomic sequences were generated from positive samples collected, along with epidemiological metadata, at a walk-up, rapid testing site in the Mission District of San Francisco, California during 22 November to 1 December, 2020, and 10-29 January 2021. Secondary household attack rates and mean sample viral load were estimated and compared across observed variants. RESULTS: A total of 12 124 tests were performed yielding 1099 positives. From these, 928 high-quality genomes were generated. Certain viral lineages bearing spike mutations, defined in part by L452R, S13I, and W152C, comprised 54.4% of the total sequences from January, compared to 15.7% in November. Household contacts exposed to the "California" or "West Coast" variants (B.1.427 and B.1.429) were at higher risk of infection compared to household contacts exposed to lineages lacking these variants (0.36 vs 0.29, risk ratio [RR] = 1.28; 95% confidence interval [CI]: 1.00-1.64). The reproductive number was estimated to be modestly higher than other lineages spreading in California during the second half of 2020. Viral loads were similar among persons infected with West Coast versus non-West Coast strains, as was the proportion of individuals with symptoms (60.9% vs 64.3%). CONCLUSIONS: The increase in prevalence, relative household attack rates, and reproductive number are consistent with a modest transmissibility increase of the West Coast variants. Summary: We observed a growing prevalence and modestly elevated attack rate for "West Coast" severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants in a community testing setting in San Francisco during January 2021, suggesting its modestly higher transmissibility.

Estimation of secondary household attack rates for emergent SARS-CoV-2 variants detected by genomic surveillance at a community-based testing site in San Francisco.

BACKGROUND: Sequencing of the SARS-CoV-2 viral genome from patient samples is an important epidemiological tool for monitoring and responding to the pandemic, including the emergence of new mutations in specific communities. METHODS: SARS-CoV-2 genomic sequences were generated from positive samples collected, along with epidemiological metadata, at a walk-up, rapid testing site in the Mission District of San Francisco, California during November 22-December 2, 2020 and January 10-29, 2021. Secondary household attack rates and mean sample viral load were estimated and compared across observed variants. RESULTS: A total of 12,124 tests were performed yielding 1,099 positives. From these, 811 high quality genomes were generated. Certain viral lineages bearing spike mutations, defined in part by L452R, S13I, and W152C, comprised 54.9% of the total sequences from January, compared to 15.7% in November. Household contacts exposed to "West Coast" variants were at higher risk of infection compared to household contacts exposed to lineages lacking these variants (0.357 vs 0.294, RR=1.29; 95% CI:1.01-1.64). The reproductive number was estimated to be modestly higher than other lineages spreading in California during the second half of 2020. Viral loads were similar among persons infected with West Coast versus non-West Coast strains, as was the proportion of individuals with symptoms (60.9% vs 64.1%). CONCLUSIONS: The increase in prevalence, relative household attack rates, and reproductive number are consistent with a modest transmissibility increase of the West Coast variants; however, additional laboratory and epidemiological studies are required to better understand differences between these variants. SUMMARY: We observed a growing prevalence and elevated attack rate for "West Coast" SARS-CoV-2 variants in a community testing setting in San Francisco during January 2021, suggesting its modestly higher transmissibility.

The antimalarial efficacy and mechanism of resistance of the novel chemotype DDD01034957.

New antimalarial therapeutics are needed to ensure that malaria cases continue to be driven down, as both emerging parasite resistance to frontline chemotherapies and mosquito resistance to current insecticides threaten control programmes. Plasmodium, the apicomplexan parasite responsible for malaria, causes disease pathology through repeated cycles of invasion and replication within host erythrocytes (the asexual cycle). Antimalarial drugs primarily target this cycle, seeking to reduce parasite burden within the host as fast as possible and to supress recrudescence for as long as possible. Intense phenotypic drug screening efforts have identified a number of promising new antimalarial molecules. Particularly important is the identification of compounds with new modes of action within the parasite to combat existing drug resistance and suitable for formulation of efficacious combination therapies. Here we detail the antimalarial properties of DDD01034957-a novel antimalarial molecule which is fast-acting and potent against drug resistant strains in vitro, shows activity in vivo, and possesses a resistance mechanism linked to the membrane transporter PfABCI3. These data support further medicinal chemistry lead-optimization of DDD01034957 as a novel antimalarial chemical class and provide new insights to further reduce in vivo metabolic clearance.

Evaluation of whole genome amplification and bioinformatic methods for the characterization of Leishmania genomes at a single cell level.

Here, we report a pilot study paving the way for further single cell genomics studies in Leishmania. First, the performances of two commercially available kits for Whole Genome Amplification (WGA), PicoPLEX and RepliG were compared on small amounts of Leishmania donovani DNA, testing their ability to preserve specific genetic variations, including aneuploidy levels and SNPs. We show here that the choice of WGA method should be determined by the planned downstream genetic analysis, PicoPLEX and RepliG performing better for aneuploidy and SNP calling, respectively. This comparison allowed us to evaluate and optimize corresponding bio-informatic methods. As PicoPLEX was shown to be the preferred method for studying single cell aneuploidy, this method was applied in a second step, on single cells of L. braziliensis, which were sorted by fluorescence activated cell sorting (FACS). Even sequencing depth was achieved in 28 single cells, allowing accurate somy estimation. A dominant karyotype with three aneuploid chromosomes was observed in 25 cells, while two different minor karyotypes were observed in the other cells. Our method thus allowed the detection of aneuploidy mosaicism, and provides a solid basis which can be further refined to concur with higher-throughput single cell genomic methods.

Complete Genome Sequence of a Novel Coronavirus (SARS-CoV-2) Isolate from Bangladesh.

The complete genome sequence of a novel coronavirus (severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2]) isolate obtained from a nasopharyngeal swab from a patient with COVID-19 in Bangladesh is reported.

Inhibition of Resistance-Refractory P. falciparum Kinase PKG Delivers Prophylactic, Blood Stage, and Transmission-Blocking Antiplasmodial Activity.

The search for antimalarial chemotypes with modes of action unrelated to existing drugs has intensified with the recent failure of first-line therapies across Southeast Asia. Here, we show that the trisubstituted imidazole MMV030084 potently inhibits hepatocyte invasion by Plasmodium sporozoites, merozoite egress from asexual blood stage schizonts, and male gamete exflagellation. Metabolomic, phosphoproteomic, and chemoproteomic studies, validated with conditional knockdown parasites, molecular docking, and recombinant kinase assays, identified cGMP-dependent protein kinase (PKG) as the primary target of MMV030084. PKG is known to play essential roles in Plasmodium invasion of and egress from host cells, matching MMV030084's activity profile. Resistance selections and gene editing identified tyrosine kinase-like protein 3 as a low-level resistance mediator for PKG inhibitors, while PKG itself never mutated under pressure. These studies highlight PKG as a resistance-refractory antimalarial target throughout the Plasmodium life cycle and promote MMV030084 as a promising Plasmodium PKG-targeting chemotype.

Probing the Open Global Health Chemical Diversity Library for Multistage-Active Starting Points for Next-Generation Antimalarials.

Most phenotypic screens aiming to discover new antimalarial chemotypes begin with low cost, high-throughput tests against the asexual blood stage (ABS) of the malaria parasite life cycle. Compounds active against the ABS are then sequentially tested in more difficult assays that predict whether a compound has other beneficial attributes. Although applying this strategy to new chemical libraries may yield new leads, repeated iterations may lead to diminishing returns and the rediscovery of chemotypes hitting well-known targets. Here, we adopted a different strategy to find starting points, testing ∼70,000 open source small molecules from the Global Health Chemical Diversity Library for activity against the liver stage, mature sexual stage, and asexual blood stage malaria parasites in parallel. In addition, instead of using an asexual assay that measures accumulated parasite DNA in the presence of compound (SYBR green), a real time luciferase-dependent parasite viability assay was used that distinguishes slow-acting (delayed death) from fast-acting compounds. Among 382 scaffolds with the activity confirmed by dose response (<10 µM), we discovered 68 novel delayed-death, 84 liver stage, and 68 stage V gametocyte inhibitors as well. Although 89% of the evaluated compounds had activity in only a single life cycle stage, we discovered six potent (half-maximal inhibitory concentration of <1 µM) multistage scaffolds, including a novel cytochrome bc1 chemotype. Our data further show the luciferase-based assays have higher sensitivity. Chemoinformatic analysis of positive and negative compounds identified scaffold families with a strong enrichment for activity against specific or multiple stages.

Combining Stage Specificity and Metabolomic Profiling to Advance Antimalarial Drug Discovery.

We report detailed susceptibility profiling of asexual blood stages of the malaria parasite Plasmodium falciparum to clinical and experimental antimalarials, combined with metabolomic fingerprinting. Results revealed a variety of stage-specific and metabolic profiles that differentiated the modes of action of clinical antimalarials including chloroquine, piperaquine, lumefantrine, and mefloquine, and identified late trophozoite-specific peak activity and stage-specific biphasic dose-responses for the mitochondrial inhibitors DSM265 and atovaquone. We also identified experimental antimalarials hitting previously unexplored druggable pathways as reflected by their unique stage specificity and/or metabolic profiles. These included several ring-active compounds, ones affecting hemoglobin catabolism through distinct pathways, and mitochondrial inhibitors with lower propensities for resistance than either DSM265 or atovaquone. This approach, also applicable to other microbes that undergo multiple differentiation steps, provides an effective tool to prioritize compounds for further development within the context of combination therapies.

Genomes of Leishmania parasites directly sequenced from patients with visceral leishmaniasis in the Indian subcontinent.

Whole genome sequencing (WGS) is increasingly used for molecular diagnosis and epidemiology of infectious diseases. Current Leishmania genomic studies rely on DNA extracted from cultured parasites, which might introduce sampling and biological biases into the subsequent analyses. Up to now, direct analysis of Leishmania genome in clinical samples is hampered by high levels of human DNA and large variation in parasite load in clinical samples. Here, we present a method, based on target enrichment of Leishmania donovani DNA with Agilent SureSelect technology, that allows the analysis of Leishmania genomes directly in clinical samples. We validated our protocol with a set of artificially mixed samples, followed by the analysis of 63 clinical samples (bone marrow or spleen aspirates) from visceral leishmaniasis patients in Nepal. We were able to identify genotypes using a set of diagnostic SNPs in almost all of these samples (97%) and access comprehensive genome-wide information in most (83%). This allowed us to perform phylogenomic analysis, assess chromosome copy number and identify large copy number variants (CNVs). Pairwise comparisons between the parasite genomes in clinical samples and derived in vitro cultured promastigotes showed a lower aneuploidy in amastigotes as well as genomic differences, suggesting polyclonal infections in patients. Altogether our results underline the need for sequencing parasite genomes directly in the host samples.

Open-source discovery of chemical leads for next-generation chemoprotective antimalarials.

To discover leads for next-generation chemoprotective antimalarial drugs, we tested more than 500,000 compounds for their ability to inhibit liver-stage development of luciferase-expressing Plasmodium spp. parasites (681 compounds showed a half-maximal inhibitory concentration of less than 1 micromolar). Cluster analysis identified potent and previously unreported scaffold families as well as other series previously associated with chemoprophylaxis. Further testing through multiple phenotypic assays that predict stage-specific and multispecies antimalarial activity distinguished compound classes that are likely to provide symptomatic relief by reducing asexual blood-stage parasitemia from those which are likely to only prevent malaria. Target identification by using functional assays, in vitro evolution, or metabolic profiling revealed 58 mitochondrial inhibitors but also many chemotypes possibly with previously unidentified mechanisms of action.

Defining the Determinants of Specificity of Plasmodium Proteasome Inhibitors.

The Plasmodium proteasome is an emerging antimalarial target due to its essential role in all the major life cycle stages of the parasite and its contribution to the establishment of resistance to artemisinin (ART)-based therapies. However, because of a similarly essential role for the host proteasome, the key property of any antiproteasome therapeutic is selectivity. Several parasite-specific proteasome inhibitors have recently been reported, however, their selectivity must be improved to enable clinical development. Here we describe screening of diverse libraries of non-natural synthetic fluorogenic substrates to identify determinants at multiple positions on the substrate that produce enhanced selectivity. We find that selection of an optimal electrophilic "warhead" is essential to enable high selectivity that is driven by the peptide binding elements on the inhibitor. We also find that host cell toxicity is dictated by the extent of coinhibition of the human ß2 and ß5 subunits. Using this information, we identify compounds with over 3 orders of magnitude selectivity for the parasite enzyme. Optimization of the pharmacological properties resulted in molecules that retained high potency and selectivity, were soluble, sufficiently metabolically stable and orally bioavailable. These molecules are highly synergistic with ART and can clear parasites in a mouse model of infection, making them promising leads as antimalarial drugs.

UCT943, a Next-Generation Plasmodium falciparum PI4K Inhibitor Preclinical Candidate for the Treatment of Malaria.

The 2-aminopyridine MMV048 was the first drug candidate inhibiting Plasmodium phosphatidylinositol 4-kinase (PI4K), a novel drug target for malaria, to enter clinical development. In an effort to identify the next generation of PI4K inhibitors, the series was optimized to improve properties such as solubility and antiplasmodial potency across the parasite life cycle, leading to the 2-aminopyrazine UCT943. The compound displayed higher asexual blood stage, transmission-blocking, and liver stage activities than MMV048 and was more potent against resistant Plasmodium falciparum and Plasmodium vivax clinical isolates. Excellent in vitro antiplasmodial activity translated into high efficacy in Plasmodium berghei and humanized P. falciparum NOD-scid IL-2Rγ null mouse models. The high passive permeability and high aqueous solubility of UCT943, combined with low to moderate in vivo intrinsic clearance, resulted in sustained exposure and high bioavailability in preclinical species. In addition, the predicted human dose for a curative single administration using monkey and dog pharmacokinetics was low, ranging from 50 to 80 mg. As a next-generation Plasmodium PI4K inhibitor, UCT943, based on the combined preclinical data, has the potential to form part of a single-exposure radical cure and prophylaxis (SERCaP) to treat, prevent, and block the transmission of malaria.